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We provide a straightforward and numerically efficient procedure to perform local density approx- 
imation + Hubbard I (LDA+HIA) calculations, including self-consistency over the charge density, 
within the full potential linearized augmented plane wave (FP-LAPW) method. This implemen- 
tation is all-electron, includes spin-orbit interaction, and makes no shape approximations for the 
charge density. The method is applied to calculate selected heavy actinides in the paramagnetic 
phase. The electronic structure and spectral properties of Am and Cm metals obtained are in 
agreement with previous dynamical mean-field theory (LDA+DMFT) calculations and with avail- 
able experimental data. We point out that the charge density self-consistent LDA+HIA calculations 
predict the / charge on Bk to exceed the atomic integer / 8 value by 0.22. 

PACS numbers: 71.27.+a,79.60.-i 



I. INTRODUCTION 

It is known that conventional band theory — 
local density approximation (LDA) and its semi- 
local extension, generalized gradient approximation 
(GGA) — gives poor results for actinides. Since the 
LDA/GGA results are qualitatively incorrect al- 
ready at the level of ground state properties, like 
the equilibrium volume and magnetization, the elec- 
tronic structure theory of actinides requires that 
electron-electron correlations arc included beyond 
those given by conventional LDA/GGA. Lately, sev- 
eral correlated band theory approaches have been 
put forward: LDA+Hubbard U (LDA+U) 0,0, the 
hybrid functional (HYF) approach Q or the self- 
interaction-corrected local spin-density (SIC-LSD) 
Q ■ Each of them achieved an improvement of some 
particular aspects of the electronic structure of ac- 
tinides. 

None of these correlated band theories has been 
capable to correctly describe spectral properties of 
actinides. Recently, the excitations in Pu and Am 
were extensively studied with the aid of a combi- 
nation of the LDA and the dynamical mean-field 
theory (LDA+DMFT) @, i, 0, i, H , that success- 
fully explains the experimentally observed multi- 
peak structure in Pu valence-band photoclcctron 
spectra (PES). In spite of obvious progress in the 
LDA+DMFT theory, it has been mostly focused on 
calculations of excitations and implemented on the 
basis of a tight-binding Hamiltonian built from the 
LDA, without self-consistency over the charge den- 



sity. 

In this paper we present a simple and numerically 
efficient procedure to combine the LDA+Hubbard I 
approximation (HIA), including self-consistency 
over the charge density, with the full potential lin- 
earized augmented plane wave (FP-LAPW) method 
0. The FP-LAPW method makes no shape ap- 
proximation for the charge density and is considered 
to be state-of-the-art in accuracy. We apply our im- 
plementation to the electronic structure and spectro- 
scopic properties of heavy actinides: Am, Cm and 
Bk. 

There is a revival of interest to the electronic and 
spectroscopic properties of heavy actinides [Tl| . Su- 
perconducting temperature of Am shows complex 
and unconventional dependence on lattice structure 
transformations ■ On the basis of standard band 
structure calculations it was proposed that curium 
is one of a few elements that has its lattice struc- 
ture stabilized by magnetism [H[. The spectro- 
scopic studies [l4[ suggested that 5/ states of Cm are 
shifted towards the LS coupling limit, unlike most 
actinide elements where the intermediate coupling 
prevails. 

The paper is organized as follows. For the sake 
of completeness, in Sec. II we recall the basic 
equations of the LDA+DMFT in a formulation of 
Ref. [l5|. Then we describe charge density self- 
consistent LDA+HIA approximation implemented 
in FP-LAPW method. In Sec. Ill we present the re- 
sults of the charge density self-consistent LDA+HIA 
calculations for 5/ Am, Cm, and Bk elemental met- 
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als in the paramagnetic state. These results are com- 
pared with previous work and additional features are 
pointed out. 



II. METHODOLOGY 

We start with the multi-band Hubbard Hamilto- 
nian [H H = H° + H int , where 



1° ■ c f c 



h° = E E H l 

i,j 71.72 

= EE<,T,( k X( k )s(k) (i) 

k 71,72 

is the one-particle Hamiltonian found from ab initio 
electronic structure calculations of a periodic crystal, 
including the spin-orbit coupling (SOC). The indices 
i, j label lattice sites, 7 = (Ima) denote spinorbitals 
{</> 7 }, and k is a k- vector from the first Brillouin 
zone. It is assumed that the electron-electron cor- 
relations between s, p, and d electrons are well de- 
scribed within the density functional theory, while 
the correlations between the / electrons have to be 
considered separately by introducing the interaction 
Hamiltonian 



H°(k) and a guess for the local selfenergy 
the local Green function is calculated by integrat- 
ing G(k, z), Eq. ([3]), over the Brillouin zone. Sub- 
sequently bath Green function (the so-called Weiss 
field) is calculated [Hj], and used to solve the SIAM. 
New local is evaluated, which is inserted back 
into Eq. ([3]). In addition, the charge density needed 
to construct the single-particle Hamiltonian Ho(k) 
in Eq. ([3J has to be calculated self-consistently from 
the local Green function. 



A. Hubbard-I approximation for E(z) 

We make use of the multiorbital HIA, which is 
suitable for incorporating the multiplet transitions 
into the electronic structure, as it is explicitly based 
on the exact diagonalization of an isolated atomic- 
like shell. Further, we restrict our formulation to the 
paramagnetic phase. In HIA, only the site-diagonal 
terms from one-particle Hamiltonian Eq. |T]) are re- 
tained [l5j . and the on-site atomic-like Hamiltonian 
including SOC is constructed, see also Ref. Q, 
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The operator V ee represents an effective on-site 
Coulomb interaction |l5[ expressed in terms of the 
Slater integrals and the spherical harmonics \lm). 

In what follows we use a local approximation for 
the one-particle selfenergy S(k, z) which contains 
the electron-electron correlations, i. e., we assume 
that the selfenergy is site-diagonal and therefore 
independent of k. The corresponding one-particle 
Green function reads 



G(k, z) = (z + ij, - H"{k) - 



(3) 



where z is a (complex) energy measured with respect 
to the chemical potential \x. The interaction term, 
Eq. @, acts only in the subspace of /-states. Con- 
sequently, the selfenergy S(z) is nonzero only in the 
subspace of the /-states. 

The self-consistent procedure to solve the peri- 
odic lattice problem in the DMFT approximation is 
now formulated in the usual way making use of the 
"impurity" method of Ref. p]|. The DMFT self- 
consistency condition is achieved by equating the 
local Green function in a solid to the Green func- 
tion of a single-impurity Anderson model (SIAM) 
that describes an isolated multiorbital impurity sur- 
rounded by a bath of uncorrelated delocalized elec- 
trons. Starting with the single-particle Hamiltonian 



+ - ^2 {mim 2 \V ee \m 3 m4) 

mi . . .7714 
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where £ is the SOC parameter. Note that the crystal 
field terms are not included in Eq. (U|), and will be 
treated on the LDA level that is sufficient for our 
applications. 

Consecutively, exact diagonalization, H &Xl \v) = 
E v \v). is performed in order to obtain all eigenvalues 
E v and eigenvectors | v) that are used to calculate the 
atomic Green function 
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Here (3 is the inverse temperature, Z is the parti- 
tion function, and N v is the number of particles in 
the state \v). These N v are eigenvalues of the parti- 
cle number operator that commutes with the atomic 
Hamiltonian Eq. (H}. Parameter /i# plays a role of 
the HIA chemical potential. Actual choice of hh will 
be discussed later. Finally, the atomic self-energy is 
evaluated as 
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This Tih(z) contains all local spin-orbit and 
Coulomb correlation effects. 



B. Self-consistency over charge density: 
Local Density Matrix Approximation 

Instead of solving Eq. ([3]) directly, we look for an 
approximate solution including charge density self- 
consistency in a way which is similar to the well 
known rotationally invariant LDA+U method (Id1 |. 

We start with calculating the HIA £(z), Eq. ©, 
for given [in. In our applications, this starting hh 
corresponds to the nominal atomic /-shell occupa- 
tion nf. The initial solution for lattice electrons is 
represented by the LDA Green function matrix in 
the local basis {</> 7 }, 
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Note that the SOC is included in the LDA Hamilto- 
nian Hloa(^)- The local impurity Green function 
is calculated combining and Glda(z), 



G(z) 
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Glda(z) 
- ( AeS. 



-l 

7172 



7172 



(8) 



where Ae is chosen to keep the given number of /- 
electrons n / , and serves as an analogon of the differ- 
ence between the impurity and the lattice chemical 
potentials fl7j |. 

With the aid of G(z) from Eq. ((8]), the occupation 
matrix n 7l72 = — 7r _1 lm J F dz [G(z)] 7l72 is evalu- 
ated , and used to construct the effective "LDA+U 
potential" V v = E 7l72 I K W 72 (</>72 1, where 



^7172 = £ ({ 727 |F ee | 7l7 ') - < 7 27l^ ee |7'7i))« 
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(9) 



In what follows, we have adopted the fully-localized 
(or atomic-like) limit (FLL) prescription of Solovyev 
et al. [2(| for the double-counting term Vdc = U(nf — 
1/2) - J(n f - l)/2. 

The set of Kohn-Sham-like equations is solved self- 
consistently over the charge density p(r) 



V 2 + V LDA (r) +Vu + £(1 • s)J $,(r) = e^r) , 

OCC 

p(r)=2*J(r)* i (r) > (10) 



where the effective potential is the sum of the 
standard LDA potential Vlda(t) and the on-site 
electron-electron interaction potential Vjj- Solving 
Eqs. (JTUJ) is similar to solving Eq. ([3]) in a sense that 
the selfenergy matrix £(z) from Eq. (J6]) is substi- 
tuted by the energy independent potential matrix 
defined in Eq. ©. After the self-consistency over 
the charge density is achieved, the LDA+U Green 
function matrix Gu{z) in the local basis {</> 7 } is cal- 
culated from Eq. ((7|), substituting Hlda by LDA+U 
Hamiltonian. Finally, new uncorrclatcd Green func- 
tion 



G 



lda{ 



(z)= G u 1 (z) + V u (z) 



(11) 



is evaluated. The self-consistency loop is closed by 
inserting this new Glda(z) into the matrix equa- 
tion ||5J). In addition, an updated selfenergy E(z) 
is calculated with the aid of Eqs. (|5l6p . where the 
new value of /j,h is set equal to the double-counting 
potential Vdc that corresponds to n/ obtained from 
the LDA+U Green function. 

The condition = Vd c is essential and can be 
justified as follows. The double-counting term Vdc 
accounts approximately for the electron-electron in- 
teraction energy Ef^ DA already included in the LDA. 
Namely, Vdc is a derivative of this mean energy con- 
tribution with respect to the /-shell occupation n/, 
Vdc = 9Ef^ DA /dnf. Indeed, it represents a mean- 
field value of the chemical potential fin that controls 
the number of / electrons. 

The FLL choice of the double counting Vd c 

is not unique and other prescriptions, for instance 
the so-called "around-mean-field" Vdc pH . [HJ , can 
be used. Up to date, there is no precise solution for 
the double counting in the conventional LDA/GGA 
as it does not have a diagrammatic representation 
that would provide explicit identification of the cor- 
responding many-body interaction terms. There- 
fore, "physical" arguments prevail in the choice of 
Vdc- Since we will be dealing with heavy actinides 
with well localized /-manifolds, it is reasonable to 
use the FLL double counting that is assumed to per- 
form better for the case of /-occupation close to in- 
teger. 

We will refer to our procedure as the "local 
density matrix approximation" (LDMA), since full 
convergence for Glda, p(r) and hh is achieved 
when the local occupation matrix n 7l72 is con- 
verged. We would like to emphasize that the self- 
consistency condition of equating the occupation 
matrix obtained from the local impurity Green func- 
tion Eq. ([5]) to the local occupation matrix in solid 
(used in the LDA+U potential Eq. (jTQJ) ) is a subset 
of general DMFT condition that the SIAM Green 
function is equal to the local Green function in a 
solid H|. 

What makes our approach different from the con- 
ventional LDA+HIA given by Eq. ^ and from sim- 
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ilar basis set extension method of Ref. [24|, is that 
we interchange the "inner" DMFT self-consistency 
loop over the bath Green function Glda, Eq. (fTT|) . 
and the "outer" self-consistency loop over the charge 
density p{v), Eq. (fT0|) . 

Up to now, our considerations did not depend on 
the choice of the basis set. The method becomes ba- 
sis dependent, when a projector for the Bloch state 
3>i(r) solution of Eq. (fT0|) on the local basis {</> 7 } 
is specified. The FP-LAPW method uses a basis 
set of plane waves that are matched onto a linear 
combination of all radial solutions (and their energy 
derivative) inside a sphere centered on each atom. 
In this case, we make use of the projector technique 
which is described in detail in Ref. It is im- 

portant to mention that due to the full potential 
character care should be taken [l8| to exclude the 
double-counting of the /-states non-spherical contri- 
butions to the LDA and LDA+U parts of potential 
in Eq. (TO]). 



III. RESULTS 

As representative systems to illustrate the LDMA 
numerical procedure, we select heavy actinides — 
Am, Cm and Bk. For all of them, the HIA is ex- 
pected to provide a reasonable approximation for 
the selfenergy. We focus on comparison between 
the theory and available experimental results for 
valence-band photoelectron spectra (PES) as well 
as X-ray absorption (XAS) and electron energy-loss 
(EELS) spectroscopies. This comparison is often 
taken as important criterion of truthfulness of elec- 
tronic structure calculations. 

Experimental valence-band PES spectra will be 
compared with valence spectral densities resulting 
from the self-consistent LDMA. For the XAS and 
EELS experiments, we will compare the branching 
ratio B as well as strength of the spin-orbit coupling 
w 110 for core-to-valence 4d-5/ transition, 

110 7/2 4 5/2 

w n f ~ 3 n / ' 

where A represents a small correction term [ill ] . 

The actinides were calculated assuming paramag- 
netic state with /cc-crystal structure and the exper- 
imental volume per atom. The parameters of the 
local Hamiltonian, Eq. (j4|), were chosen as follows: 
U = F = 4.5 eV, F<2, F4 and F 6 were taken from 
Table III of Ref. [ll(, and values of the SOC param- 
eter £ were extracted from LDA calculations = 
0.35 eV, icm = 0.36 eV and £ Bk = 0.42 eV). The 
HIA Green function and selfenergy, Eqs. ((5]) and (J6j> , 
were calculated along the real axis z = Rez + iS 



with S = 0.1 eV. In the process, values 10 eV -1 and 
100 eV -1 were used for the inverse temperature (3. 
For self-consistency, 108 special k-points [25| in the 
irreducible l/8th part of the BZ were used. The 
same sphere radius Rmt = 3.1 a.u. was used for 
all actinides, and Rmt x K max = 10.70 determined 
the basis set size. The /-manifold occupation nt is 
varied in the calculations until the convergence bet- 
ter than 0.01 for rif and 0.001 for all components 
of the on-site occupation matrix n 7l72 is achieved. 
The charge density is fully converged to better than 
10~ 5 c/a.u. 3 at each iteration. 

We plot in Fig. [1] total and /-projected spectral 
densities resulting from self-consistent LDMA calcu- 
lations (i. c., from converged Eq. (JSJ)). In the case of 
Am, wc obtain very good agreement with previous 
LDA+DMFT calculations [2I HH as we ll as with 
our own non-self-consistent LDA+HIA calculations 
Q for spectral peak positions in occupied and unoc- 
cupied parts of the spectrum. Comparison with the 
PES experimental data [27| is very good. No PES 
and BIS measurements exist for Cm and Bk metals. 
The results of present Cm calculations agree reason- 
ably well with the results of recent DMFT study @ 
supporting validity of the LDMA. We found practi- 
cally no changes in the densities of states when f3 was 
decreased from 10 cV^ 1 to 100 eV -1 . The results 
turn out to be almost insensitive to the choice of /?, 
since it enters explicitly only Eq. §5§ and has prac- 
tically no influence on the chemical potential fin- 

For Am and Cm, the self-consistent value of n/ is 
very close to the atomic integer value (see Tab. [XJ) in 
agreement with LDA+DMFT results @, Hi]. For 
Bk, deviation of n/ from the nominal atomic / 8 
is somewhat bigger (Tab. suggesting a possibil- 
ity of mixed- valence character in some of heavy ac- 
tinides. In fact, Svane et al. Q have already sug- 
gested mixed-valence states in Am, Cm and Bk 
on the basis of SIC-LSD calculations (see Tab. I 
of Ref. [3]) that split the /-electrons into localized 
manifold with fixed valence and an itinerant part. 
Present LDMA as well as previous LDA+DMFT 
d, [24| calculations show that the tendency to mixed- 
valence in heavy actinides is substantially overesti- 
mated by the SIC-LSD theory. 

Now we turn to comparison with XAS and EELS 
experiments [l4|. In these experiments, the inten- 
sities 7 5/2 (4d 5/2 -> 5/5/2,7/2) and I 3/2 (4d 3/2 -> 
5/5/2) of the X-ray absorption lines are measured 
and the branching ratio B = 15/2/(^3/2 + I5/2) is 
obtained. Note that B is the only quantity which di- 
rectly follows from the experiments. To extract the 
SOC strength w 110 , the atomic sum rules are used in 
conjunction with the atomic calculations [Tlj . In or- 
der to compare with the experiment, we obtain 715/2 
and 717/2 from the local occupation matrix n 7l72 and 
make use of Eq. (Tj"2"j) to obtain B and w 110 . We do 
not take into account the small correction factor A 
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The LDMA results for Am are shown in Table [J 
in comparison with the experimental data [28| and 
the results of atomic intermediate-coupling (IC) cal- 
culations [Hi- The LDMA calculated 715/2, "7/2, 
branching ratio B, and spin-orbit coupling strength 
are close to atomic IC and experimentally derived 
values. Once again, present calculations confirm lo- 
calized nature of solid state Am /-manifold close to 
the atomic / 6 configuration. 

The LDMA results for Cm are also shown in Ta- 
ble U in comparison with the results of DMFT calcu- 
lations [|| , atomic IC calculations [2!| as well as with 
experimental data [l4[. There is a very good agree- 
ment for n 5 /2, 71.7/2, B, and w 110 /(14 — nf) between 
LDMA and atomic IC calculations. Also, the calcu- 
lated branching ratio agrees with B=0.75 obtained 
from DMFT calculations Q. Note that LDMA, 
DMFT and IC results agree with each other, and 
all slightly differ from experimentally observed B 
of 0.794 []J|. Gaining inspiration from LDA/GGA, 
Moore et al. [3] suggested that Cm 5 /-states are 
shifted towards the LS coupling limit due to en- 
hancement of the exchange interaction over the spin- 
orbit coupling. However, Shim et al. [3(| noticed 
very recently that agreement between the theory and 
experiment for B improves substantially when the 
Slater integrals [ll[ are slightly reduced to account 
for the solid state screening. 

To date, no XAS or EELS experimental data exist 
for Bk metal. The calculated "5/2, 17/2, branching 
ratio B, and spin-orbit coupling strength arc listed 
in Table [T] together with the atomic IC / 8 calcula- 
tions. The main difference between the solid state 
and the atomic /-manifolds is due to an increase in 
occupation of 77,7/2 states. Nevertheless, the values 
of B and w 110 per hole are practically the same. 
The measurements of the branching ratio are often 
used to obtain the experimental value of nf. Our 
results illustrate that the knowledge of the B-ratio 
alone is not sufficient for precise determination of 
the /-manifold occupation. 

Now we turn to an estimate of the effective local 
magnetic moment (i e ff i n the paramagnetic phase. 
Importance of magnetism in Cm metal was empha- 
sized recently in the context of its phase stability 
[l3j |. The temperature independent magnetic sus- 
ceptibility is found for Am [3l[ that is consistent 
with zero /i e //. The magnetic susceptibility mea- 
surements in the paramagnetic phase yield effective 
magnetic moment of ~ 8 (is for Cm and ~ 9.8 (Lb 
for Bk [H|. 

We can estimate semi-quantitativcly the effective 
local moment making use of the atomic Hamiltonian, 
Eq. (|4]), and the chemical potential (in = Vdc that 
is self-consistently determined in the LDMA calcu- 
lations. The expectation values of total, spin and 
orbital moment operators, J, S and L, are calcu- 



Am 


n f 


5/2 
n f 


7/2 
n,f 


B 


w 110 /n h 


LDMA (/3= 10 eV _1 ) 


5.95 


5.11 


0.83 


0.897 


-0.743 


LDMA ((3= 100 eV _1 ) 


5.95 


5.16 


0.79 


0.902 


-0.756 


atomic IC fill 


6 


5.28 


0.72 


0.916 


-0.79 


Exp. [14] 


6 


5.38 


0.62 


0.930 


-0.825 


Cm 


n f 


5/2 
rif 


772 
n f 


B 


w 110 /n h 


LDMA ((3= 10 eV ) 


7.07 


4.04 


3.03 


0.736 


-0.340 


LDMA ((3= 100 eV" 1 ) 


7.07 


4.04 


3.03 


0.737 


-0.341 


DMFT [6] 


7.0 


N/A N/A 


0.75 


N/A 


atomic IC fill 


7 


4.10 


2.90 


0.75 


-0.37 


Exp. [14] 


7 


4.41 


2.59 


0.794 


-0.485 


Bk 


71/ 


5/2 
n f 


7/2 
n f 


B 


w lw /n h 


LDMA ((3= 10 eV _i ) 


8.22 


5.01 


3.21 


0.840 


-0.591 


LDMA {(3= 100 eV" 1 ) 


8.22 


5.01 


3.21 


0.840 


-0.601 


atomic IC [11] 


8 


5.00 


3.00 


0.84 


-0.61 



TABLE I: Branching ratio B and spin-orbit coupling 



u /rih, where n/, = (14 — n/), for 
Note that "experimental" values of 



strength per hole 
Am, Cm and Bk. 

n°f 2 and n^ 2 are not measured, but derived from sum 
rule Eq. (|12|l assuming integer atomic occupation nf. 



lated as grand-canonical averages, 



<X'> = -Tr 



X. 2 eM-l3[H at -^HN}) 



J,S,L 



(13) 



Further, spin S, orbital L and total J moment 
"quantum numbers" are found using (X 2 ) = X(X + 
1) for X = S, L, J. Subsequently, the effective mag- 
netic moment /i e // = gj\J J{J + 1) is evaluated, 
where the g- factor gj = (2S + L)/J is used. 

For Am, we obtain S = L = 2.33 and J = for 
(3 = 100 cV _1 in Eq. (p~3|) . Decreasing the value of 
(3 to 10 eV -1 yields a small difference in S and L 
values, and gives a non-zero value of J = 0.099. It 
means that the thermal population of the multiplcts 
excited over the non-magnetic / 6 ground state starts 
to produce non-negligible contribution in Eq. (|13[) . 

For Cm, S = 3.30, L = 0.40, and J = 3.50 are cal- 
culated from Eq. ^ ior (3 = 100 eV -1 . Decrease of 
(3 to 10 eV -1 produces practically no difference in S, 
-L, and J values. The corresponding local magnetic 
moment (J- e ff = 7.94 \xb agrees well with atomic IC 
value and experimental data [3l| (see Tab. HI)) . 

For Bk, Eq. (13]) yields S = 2.71, L = 0.40, and 
J = 6.00 for (3 = 100 eV" 1 , as well as for (3 = 
10 eV -1 . The effective magnetic moment /i e // = 
9.8 (Lb agrees well with the atomic / 8 IC value and 
experimental data [3l[ shown in Table ITT1 

Our calculations, which are not bound by any par- 
ticular atomic coupling scheme, illustrate once again 
that IC scheme is suitable for heavy actinides. Also, 
a good agreement of estimated fi e f f with experimen- 
tal data is somewhat surprising. 
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fJ-eff (MS) 


Am 


Cm 


Bk 


LDMA 





7.94 


9.54 


IC [31] 





7.6 


9.3 


Exp. [31] 





~ 8 


~ 9.8 



TABLE II: Effective local magnetic moment (i e ff for 
Am, Cm and Bk. The atomic IC vales of /J, e ff and ex- 
perimental data [3l[ are shown. 



IV. DISCUSSION AND CONCLUSIONS 

For a better insight, it is useful to point out that 
in the current implementation, which is based on a 
single-site approximation, Eq. ([5]) , to the solution of 
Eq. ([3]), the LDMA can be regarded as an exten- 
sion of the LDA+U. Importantly, the on-site occu- 
pation matrix n 7l72 is now evaluated in a many-body 
Hilbert space rather than in a single-particle Hilbert 
space as in the conventional LDA+U [lj|. Current 
implementation can be further extended towards a 
fully self-consistent DMFT making use of Wannicr- 
like basis set together with more sophisticated ap- 
proximation for the quantum impurity solver along 
the lines proposed in Ref. (32j. 

Our approach to the charge self-consistency is es- 
sentially different from the one proposed by Lecher- 
mann et al. [33j . The on-site occupation matrix, 
instead of the full charge density, is obtained from 
the local Green function. The corresponding orbital- 
dependent effective potential is used in Eq. (p~0|) 



to calculate a new bath Green function Glda, 
Eq. (fTTj) , instead of orbital-independent Kohn-Sham 
(LDA/GGA) potential. 

In this paper, we do not address the very impor- 
tant issues of the total energy calculation and deter- 
mination of the equilibrium lattice properties. The 
practical implementation of accurate total energy 
calculations is ongoing work that will be discussed 
in detail in the future. 

To summarize, we have presented a straightfor- 
ward and numerically efficient local density matrix 
approximation (LDMA) to perform the LDA+HIA 
calculations in the FP-LAPW basis, including self- 
consistency over the charge density. This implemen- 
tation is all-electron, incorporates spin-orbit inter- 
action, and includes no shape approximations for 
the charge density. The method works well for the 
electronic spectrum of representative actinide Am, 
Cm and Bk metals. Importantly, the method al- 
lows fully self-consistent calculations for the para- 
magnetic phase of the local moment systems with 
strong Coulomb correlations. It can be extended to 
incorporate the total energy and to treat the mag- 
netically ordered phases. 
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FIG. 1: Total-DOS and /-DOS for /co-Am, /co-Cm, and 
/cc-Bk for /?= 100 eV" 1 . 



